LOCAL-TO-GLOBAL PRINCIPLES FOR ROTOR WALK 
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f— ^ Abstract. In rotor walk on a finite directed graph, the exits from 

^SJ each vertex follow a prescribed periodic sequence. Here we consider 

w\ the case of rotor walk where a particle starts from a designated 

^ source vertex and continues until it hits a designated target set, 

.^^ at which point the walk is restarted from the source. We show 

that the sequence of successively hit targets, which is easily seen 

^^ to be eventually periodic, is in fact periodic. We show moreover 

I— I that reversing the periodic patterns of all rotor sequences causes 

^^ the periodic pattern of the hitting sequence to be reversed as well. 

r) The proofs involve a new notion of equivalence of rotor configura- 

• tions, and an extension of rotor walk incorporating time-reversed 

"^ particles. 

B 

C\l 1. Introduction 

> 

CN A rotor walk in a graph G is a walk in which the sequence of exits 

^ from each vertex is periodic. The sequence of exits from a vertex v 

t:J- is called the rotor mechanism at v. Rotor walks have been studied in 

(^ combinatorics as deterministic analogues of random walks, in computer 

O science as a means of load-balancing and territory exploration, and 

^ in statistical physics as a model of self-organized criticality. In this 

;:• paper, we explore several properties of the rotor mechanism that imply 

• ^ corresponding properties of the hitting sequence when G comes with a 

rS set of designated target vertices: 

^ • Given a (periodic) rotor mechanism at each vertex, the hitting 

sequence of the associated rotor walk is periodic (Theorem 1). 
• If every rotor mechanism is palindromic, then the hitting se- 
quence is palindromic (Theorem 3). 
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• If every rotor mechanism is m-repetitive, then the hitting se- 
quence is ?Ti-repetitive (Theorem 4). 

See below for precise definitions. Since the rotor mechanisms are 
local features of the walk — each one depends only on the exits from a 
particular vertex — while the hitting sequence is a global feature, we 
regard these theorems as local-global principles. 

Let G = (y, E) be a strongly connected finite directed graph, with 
self-loops and multiple edges permitted. For a vertex v & V , let d{y) 
denote the outdegree of v. A rotor mechanism at v is an ordering of the 
directed edges (or "arcs") emanating from w , say as e^ for 1 < z < d{y). 
Let u* denote the endpoint of the arc e^. We extend the definition of e\, 
and f * to alH G Z by taking e^ and f * to be periodic in i with period 
d{v). We often indicate the rotor mechanism at v using the notation 

V -)■ f\f^,. .. ,v'^^'"\. . . (period d{v)). 

Given a rotor mechanism at each vertex v, a rotor walk on G is a 
finite or infinite sequence of vertices Xo,Xi,X2,... in which the i-th 
occurrence of v is followed immediately by an occurrence of f*. For 
example, if the vertex set of G is {1, 2, 3, 4, 5} and the rotor mechanisms 
are 



1 


— )■ 


3,4,5,... (period 3) 


2 


— )■ 


3, . . . (period 1) 


3 


— )■ 


4, 2, . . . (period 2) 


4 


-)■ 


1, . . . (period 1) 


5 


— )■ 


1, . . . (period 1) 



then the rotor walk starting from 1 is 

1,3,4,1,4,1,5,1,3,2,3,4,1,4,1,5,1,3,2,3,4,1,... 

which is eventually periodic with period 9. Note that this sequence is 
not itself periodic (the initial 1 does not repeat) but if we isolate the 
terms equal to 4 or 5 we obtain the sequence 

4,4,5,4,4,5,... 

which is periodic with period 3. 

In general, we assume that G comes with a designated source vertex s 
that serves as the starting point of the rotor walk and a non-empty set T 
of designated target vertices such that all arcs emanating from a target 
vertex t G T go to s. In the example above, s = 1 and T = {4,5}. 
A rotor walk starts at the source vertex and always returns to the 
source vertex immediately after visiting a target vertex. We define the 
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hitting sequence as the subsequence of the rotor walk consisting of the 
terms that belong to T. In the above example, the hitting sequence is 
4 4 5 4 4 5 

It is easy to show that the hitting sequence is infinite (that is, the set 
of targets is visited infinitely often); see Lemma 6, below. As explained 
in §2, it is also easy to show that the hitting sequence is eventually 
periodic. Our first main result goes further: 

Theorem 1. The hitting sequence determined by a (periodic) rotor 
mechanism is periodic. 

As we have already seen, the rotor walk itself is typically not periodic. 
Let 7„ be the portion of the walk strictly between the n-th and (n-l-l)-st 
visits to T. In the preceding example, the sequence {7„}n>i is 

13,1,1,1323,1,1,1323,1,1,.... 

In general, this sequence is eventually periodic but is not periodic. 

A natural question is how to determine the period of the hitting 
sequence. We will see that this period divides the order of a certain 
element of the sandpile group S{G/T) of the graph G with the target 
set T collapsed to a single vertex (Lemma 20). 

Our second main result states that if we reverse the rotor mechanism 
at each vertex by replacing 

t; -^ w\ w^, . . . , v'^^''\ . . . (period d{v)) 

by 

V -^ v'^^''\v'^^''^-\ ...,v\... (period d{v)) 

for each vertex v, then the hitting sequence undergoes an analogous 
reversal; specifically, if the original hitting sequence has period D, then 
the new hitting sequence will also have period D, and for alll < i < D 
the i-th term of the new hitting sequence will equal the {D -\-l — i)-th 
term of the original hitting sequence. That is: 

Theorem 2. Reversing the periodic pattern of all rotor mechanisms 
results in reversing the periodic pattern of the hitting sequence. 

E.g., for the above example, the reversed rotor mechanism 

1 -^ 5,4,3,... (period 3) 
(period 1) 



2-^3,. 

3 ^ 2,4 

4 ^ 1,. 
5^1,. 



. . (period 2) 
(period 1) 
(period 1) 
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gives the reversed hitting sequence 5, 4, 4, 5, 4, 4, ... . 

An immediate corollary of Theorem 2 is that if the rotors are all 
palindromic (that is, if each fundamental period of each rotor reads 
the same backwards and forwards) then the same is true of the hitting 
sequence. 

Theorem 3. // all rotor mechanisms are palindromic, then the hitting 
sequence is palindromic. 

One can think of the entire collection of rotor mechanisms on G as 
a single rotor — perhaps embedded as a component of a larger system 
— whose rotor mechanism is the hitting sequence. From this perspec- 
tive, Theorems 1 and 3 are local-to-global principles asserting that if 
the sequence of exits from each vertex possesses a certain property 
(periodicity, palindromicity) , then the hitting sequence has the same 
property. We now state one more result of this type. Theorem 4. Fur- 
ther examples of local-global principles include Lemma 6, below, and 
[12, Theorem 1]. 

Call a sequence {ui}i>i m-repetitive if it consists of blocks of m 
consecutive equal terms; that is, 

"^am+l Uam+2 • • • Uam+m 

for all a > 0. 

Theorem 4. // all rotor mechanisms are m-repetitive, then the hitting 
sequence is m-repetitive. 

For example, consider the 2-repetitive rotor mechanism 

1 -^ 3,3,2,2,... (period 4) 

2 -^ 1,1,4,4,... (period 4) 

3 — )• 1, . . . (period 1) 

4 —7- 1, . . . (period 1) 

with source 1 and targets 3 and 4. The sequence of paths 7„ taken by 
the walker until it hits a target 

13, 13, 121213, 13, 124, 124, . . . (period 6) 

is not 2-repetitive, but the hitting sequence 

3,3,3,3,4,4,... (period 6) 

is 2-repetitive. 

The proof of Theorem 4 is not difficult (see §2) and uses only the 
abelian property of rotor walk (Lemma 7). The proofs of Theorems 1-3 
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make essential use of a new notion of equivalence of rotor configura- 
tions. We summarize the highlights here, referring the reader to §3.1 
for the full definitions. 

Let Vo = V — T. A rotor configuration is a map p : Vo ^ E such 
that p{v) is an arc emanating from v; the arc p{v) represents the arc 
by way of which a particle most recently exited vertex v. A particle 
configuration is a map cr : Vq — )■ N; we interpret a{v) as the number 
of particles present at vertex v. Following [11] we define an action 
{a, p) 1-^ ap of particle configurations on rotor configurations. We then 
define rotor configurations pi and p2 to be equivalent, written pi = p2, 
if there exists a particle configuration a such that api = ap2 (Lemma 10 
will show that this is an equivalence relation. In fact, pi = p2 if and 
only if o"pi = crp2 for all "sufficiently large" a, in a sense made precise 
by part (d) of Lemma 10.) We define an operation called complete 
cycle pushing which takes an arbitrary rotor configuration p as input 
and produces an acyclic rotor configuration p'^ as output. 
Theorem 5. 

(i) Each equivalence class of rotor configurations contains a unique 

acyclic configuration. 
(ii) The unique acyclic configuration equivalent to p is ep, where e 

is the recurrent identity element of the sandpile group. 
(iii) ep = p^ is the result of performing complete cycle pushing on p. 

To see the relevance of this notion of equivalence to Theorem 1, let p„ 
be the rotor configuration immediately after the rotor walk hits the tar- 
get set T for the n-th time. The sequence {p„}„>o is not periodic, but 
we will show that the sequence of equivalence classes [po], [pi], [P2], • • • 
is periodic. We then show that which target is hit by rotor walk start- 
ing at s with rotor configuration p depends only on the equivalence 
class [p]. 

In the proof of Theorem 2, a helpful trick is the use of antiparticles 
that behave like the "holes" considered in [10]: while a particle at 
vertex v first increments (progresses) the rotor at v and then moves to a 
neighbor according to the updated rotor, an antiparticle at v first moves 
to a neighbor according to the current rotor at v and then decrements 
(regresses) the rotor at v. Reversing the rotor mechanism at each vertex 
is equivalent to replacing all particles by antiparticles and vice versa. 

Related Work. Rotor walk was first studied in computer science from 
the point of view of autonomous agents patrolling a territory [18], and 
in statistical physics as a model of self-organized criticality [16]. It is 
an example of a "convergent game" of the type studied by Eriksson 
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[9] and more generally of an abelian network of communicating finite 
automata. Abelian networks were proposed by Dliar [], and their 
theory is developed in [3]. 

Rotor walk on G reflects certain features of random walk on G [ ]. 
For vertices v,w & V let d{v, w) be the number of arcs from v to w, 
and consider the Markov chain on state space V in which the transition 
probability from v to w equals d{v , w) / d{v) . The frequency Pi with 
which a particular target vertex t occurs in the hitting sequence for 
rotor walk equals the probability that the Markov chain when started 
from the source s reaches t before it reaches any other target vertex. A 
main theme of [12] and the companion article [17] is that the "global" 
discrepancy between npi and the number of times the rotor walk hits 
the target t in the first n runs is bounded — independently oi n — by 
a sum of "local discrepancies" associated with the rotors. 

A special case of the periodicity phenomenon was noted by Angel 
and Holroyd. If G is the 6-regular tree of height h and T is the set of 
leaves, it follows from the proof of Theorem 1.1 of [11] that the hitting 
sequence from the root is eventually periodic with period ^^T, and its 
fundamental period is a permutation of T. In Proposition 22 of [ ] 
Angel and Holroyd prove that for any initial setting of the rotors, the 
first #T terms of the hitting sequence are in fact a permutation of T. 



2. Abelian property, monoid action and group action 

This section collects the results from the literature that we will use. 
All of these can be found in the survey [11], and many date from 
considerably earlier; where we know of an earlier reference, we indicate 
that as well. Let d{v, w) be the number of arcs from f to w in the 
finite directed graph G. Let d{v) = ^^gy '^('^' ^) ^^ ^^^ out-degree of 
V. Let s be a designated source vertex and T be a set of designated 
target vertices. Arcs emanating from target vertices play no role in our 
argument; however, it can be helpful to imagine that for every target 
t G T we have d{t) = d{t, s) = 1 (that is, each target has just one 
outgoing arc, which points back to the source). We define Vq = V — T, 
the set of non-target vertices. We allow vertices in Vq to have arcs 
pointing to s. 

In § 1 , we defined a rotor walk on G as an infinite sequence of vertices 
Xo, Xi, X2, . . . in which the i-th occurrence of v is followed immediately 
by an occurrence of w*. The proofs make use of an alternative, "stack- 
based" picture of rotor walk, which we now describe. This viewpoint 
goes back at least to [8, 19]. 
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At each vertex f is a bi-infinite stack of cards, in which each card is 
labeled by an arc of G emanating from v. The i-th card is labeled by 
the arc Cj = (f,f*). For i > 0, the i-th card in the stack represents an 
instruction for where the particle should step upon visiting vertex v for 
the i-th. time. (When i < 0, the i-th. card in the stack never gets used, 
but it is helpful to pretend that it was used in the past before the rotor 
walk began; this point of view will play an important role in the proof 
of Theorem 3.) We also have a pointer at v that keeps track of how 
many departures from v have already occurred; this pointer moves as 
time passes. When i departures from v have occurred during the rotor 
walk thus far, we represent the state of the stack and pointer as 

[. . . , ej_2, ej-i, ej|ei+i, 6^+2, Cj+s, . . . ] 

{i is at the start of the rotor walk). Arcs Cj with j < ito the left of the 
pointer constitute the "past" of the rotor (arcs previously traversed), 
while the e/s with j > i constitute the "future" of the rotor (arcs to 
be traversed on future visits to v). The arc ej is called the retrospective 
state of the rotor. It represents the most recent arc traversed from v. 
Arc Cj+i is called the prospective state of the rotor. It represents the 
next arc to be traversed from v. When the particle next exits v (along 
arc Ci+i) the pointer moves to the right and the stack at v becomes 

[. . . , ej_2, ej_i, Cj, ej_|_i|ej+2; Cj-i-S; • • • J- 

The defining property of rotor walk is that for each vertex v, the 
sequence of labels in its stack is periodic. Initially, however, we will 
not need this assumption. We use the term stack walk to describe 
the more general situation when the stack at each vertex v may be an 
arbitrary sequence of arcs emanating from v. 

We will assume throughout that the finite directed graph G is strongly 
connected; that is, for any pair of vertices v and w there exist directed 
paths from v to w and from w to v. Note that strong connectivity 
is a global property of G. In fact, it is the only non-local ingredient 
needed for our local-global principles. The next lemma provides a sim- 
ple example of how strong connectivity parlays a local property — one 
that can be checked for each stack individually — into a corresponding 
global property of the hitting sequence. 

A sequence ai, 02, . . . whose terms belong to an alphabet A is called 
infinitive if for every a & A there are infinitely many indices i such 
that Qi = a [ ]. Thus, we say that the stack at vertex v is infinitive if 
every outgoing arc from v appears infinitely often as a label. Likewise, 
the hitting sequence is infinitive if the walk hits every target t infinitely 
often. 
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Lemma 6. // all stacks are infinitive, then the hitting sequence is in- 
finitive. 

Proof. Since G is finite, the stack walk visits at least one vertex infin- 
itely often. If the walk visits v infinitely often, then since the stack 
at V is infinitive, the walk traverses every outgoing arc from v infinitely 
often, so it visits all of the out-neighbors of v infinitely often. Since 
G is strongly connected, every vertex is reachable by a directed path 
of arcs from v, so every vertex is visited infinitely often. In particular, 
the walk hits every target infinitely often. D 

2.1. Abelian property. Suppose that several indistinguishable par- 
ticles are present on vertices of G. At each moment, one has a choice 
of which particle to move; one chooses a particle, shifts the pointer in 
the stack at the corresponding vertex, and advances that particle to a 
neighboring vertex according to the instruction on the card that the 
pointer just passed. We call this procedure a firing. 

For example, if we begin with m particles at the source vertex s, 
we can repeatedly advance one of them until it hits a target, then 
repeatedly advance another particle until it too hits a target, and so 
on, until all the particles have hit (and remain at) targets. 

The following lemma is known as the abelian property of rotor- 
routing (another name for it is the "strong convergence property," fol- 
lowing Eriksson ['']). For a proof, see [S, Theorem 4.1] or [11, Lemma 
3.9]. 

Lemma 7. Starting from particle configuration a and rotor configura- 
tion p, let vi, . . . ,Vm be a sequence of firings that results in all particles 
reaching the target set. Let N{t) be the number of particles that hit 
target t. The numbers N{t) (t G T) and the final rotor configuration 
depend only on a and p; in particular, they do not depend on the se- 
quence Vi, . . . ,Vm- 

The abelian property is all that is needed to prove Theorem 4, which 
says that if every rotor mechanism is ?7i-repetitive, then the hitting 
sequence is m-repetitive. 

Proof of Theorem 4- It suffices to show for all n that if we feed mn 
particles through the system in succession (starting them at s and 
stopping them when they hit T), then the number of particles that hit 
each target is a multiple of m; for, if we know this fact for both mn 
and m{n -\- 1), then it follows that the {mn + l)-st through {mn + m)- 
th particles must all hit the same target. By the abelian property 
(Lemma 7), if we let the mn particles walk in tandem, letting each 
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particle take its i-th step before any particle takes its (i + l)-st step, 
then since each stack is m-repetitive, the particles travel in groups of 
size m, such that the particles in each group travel the same path and 
hit the same target. D 

Note that Theorem 4 did not require the stacks to be periodic. The- 
orems 1-3 certainly do require periodic stacks, so we make this assump- 
tion for the remainder of the paper. 

2.2. Action of particle configurations on rotor configurations. 

Denote by Q the set of particle configurations 

Q = {a:Vo^N} 

and by R the set of rotor configurations 

R= {p:Vo^ E\ s{p{v)) = V for all v G Vq} 

where s(e) denotes the source of the arc e. We give Q the structure of 
a commutative monoid under pointwise addition. 

Next we recall from [11] the construction of the action 

Qx R^ R. 

Associated to each vertex w G Vq is a particle addition operator Ey 
acting on the set of rotor configurations: given a rotor configuration p, 
we define Ey{p) as the rotor configuration obtained from p by adding a 
particle at v and letting it perform rotor walk until it arrives at a target. 
Lemma 7 implies that the operators Ey commute: E^Ey, = E^jEy for 
all V, w G Vq. 

Now given a particle configuration a on G, we define 



neVb 



V J 



where the product denotes composition. Since the operators E^ com- 
mute, the order of composition is immaterial. The action of particle 
configuration a on rotor configuration p is defined by ap := Ecr{p). In 
words, ap is the rotor configuration obtained from p by placing a{v) 
particles at each vertex v and letting all particles perform rotor walk 
until they hit the target set T. By Lemma 7, the order in which the 
walks are performed has no effect on the outcome. The fact that the 
operators Ey commute ensures that we have a well-defined action, that 
is, {ai + a2)p = Cri((T2p). 
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2.3. The sandpile monoid and its action on rotor configura- 
tions. A particle configuration a is called stable if 

a{v) < d{v) - 1 for all v G Vq. 

If a is not stable, we can stabilize it by repeatedly toppling unstable 
vertices: Set (Tq = a, choose a vertex f o G Vq such that cro(uo) > divo) 
and topple it by sending one particle along each outgoing arc from vq. 
The resulting configuration ai is given by 

^ C^^ ^ / ^o(w) + d{vo, w) liw ^ Vq, 

^^ \ aQ{w) + d{vQ,w) - d{vQ) ii w = Vq. 

If 0"i is not stable, choose a vertex Vi such that cri(fi) > d{vi) and 
topple it in the same way to arrive at a new configuration o"2. Strong 
connectedness ensures that after finitely many topplings we reach a 
stable configuration, which is called the stabilization of a and denoted 
a° . The stabilization a° does not depend on the order of topplings [(i]. 
Let Q° be the set of stable particle configurations. We give Q° the 
structure of a commutative monoid with the operation 

(cri,cr2) ^ (cTi +0-2)°. 

That is, we sum the configurations pointwise, and then stabilize. By 
comparing two different toppling orders to stabilize Ui + (T2 + o"3, we 
see that ((ai + (72)° + a^Y = (cti + ((T2 + cr^)°)°, which shows that this 
operation is associative. The monoid Q° is called the sandpile monoid 
of G; its structure has been investigated in [2, 4]. 

Lemma 8. For any particle configuration a and any rotor configura- 
tion p we have 

a° p = a p. 

Hence, the action of particle configurations on rotor configurations de- 
scends to an action of the sandpile monoid 

Q° xR^R. 

Proof. [11, Lemma 3.12] We compute ap by grouping the initial rotor 
moves into "batches" each consisting of d{v) moves from a vertex v. 
The net effect of a batch of rotor moves is the same as that of a toppling 
at v. the rotor at v makes a full turn, so the rotor configuration is 
unchanged, and one particle is sent along each arc emanating from v. 
After finitely many batches, we arrive at particle configuration a° with 
rotors still configured as p. Now letting each remaining particle perform 
rotor walk until reaching the target set yields the rotor configuration 
a°p. D 
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We may express Lemma 8 as a commutative diagram 

Q^R >■ Q° X R 



R 

where the top arrow is {<J,p) ^-^ {cr°,p). 

2.4. The sandpile group and its action on spanning forests. 

We say that a stable particle configuration t E Q° is reachable from a 
particle configuration a if there exists a particle configuration r' such 
that T = {t' + (t)°. We say that r is recurrent if it is reachable from 
any a G Q. 

Note that if r is recurrent, then for any a E Q there exists Ti such 
that Ti{v) > a{v) for all v and r° = r; indeed, since r is reachable from 
a, there exists r' E Q with r = {t' + a)°, and we can take ti = t' + a. 

Denote by S{G/T) the set of recurrent particle configurations. If r is 
recurrent and a is any particle configuration, then (cr + r)° is recurrent. 
That is, the set S{G/T) is an ideal of the monoid Q° . In fact, S{G/T) 
is the minimal ideal of Q°, which shows that it is an abelian group [ ]. 
This group is called the sandpile group of G relative to the target set T. 
The set T plays the role of the sink vertex in [i i] . In the terminology of 
that paper, S{G /T) is the sandpile group of the graph G/T obtained 
by collapsing T to a single vertex. 

A rotor configuration p is acyclic if the graph (V,p(Vo)) contains no 
oriented cycles (where p(Vo) = {piy) : v G Vq}). Equivalently, the 
rotors {p{v)}vi^Va forin an oriented spanning forest of G rooted at T. 

Lemma 9. [jl j] Each addition operator Ea^ acts as a permutation on the 
set Rq of acyclic rotor configurations. Thus the action of the sandpile 
monoid Q° on rotor configurations restricts to an action 

S{G/T) xRo^Ro 

of the sandpile group S{G/T) on acyclic rotor configurations. 

A further result proved in [11] is that this group action is free and 
transitive. In other words, if p and p' are two spanning forests of 
G rooted at T, then there is a unique element of the sandpile group 
a G S{G/T) such that ap = p' . In particular, the order |5'(G/T)| of 
the sandpile group equals the number of acyclic rotor configurations 
|i?o|) which is the number of spanning forests of G rooted at T. We 
will not use these facts, however, except for a brief aside (Lemma 20) 
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where we identify the period of a sequence that arises in the proof of 
Theorem 1; there we use the freeness of the action. 

The identity element e G S{G/T) is a highly nontrivial object (see 
for instance [1 1 , Figures 4-6]) and plays a role in several of our lemmas 
below. If G has an oriented cycle, then e is distinct from the identity 
element of Q° because the latter is not recurrent. 

3. Equivalence and cycle pushing 

In this section we develop a new notion of equivalence of rotor con- 
figurations and use it to prove Theorem 1. Here and throughout the 
rest of the article, we assume that the stack at each vertex v is periodic 
with period d{v). For e = el (the ith arc in the rotor mechanism at 
v) we define e"*" = e"^^^ and e~ = e^~^, where 2 + 1 and i — 1 are to be 
interpreted modulo d{v). 

3.1. Equivalence of rotor configurations. 

Lemma 10. Let p and p' be rotor configurations on G. The following 
are equivalent: 

(a) ap = op' for some particle configuration a. 

(b) ap = ap' for all recurrent configurations a G S{G/T). 

(c) ep = ep' , where e is the identity element of S{G/T). 

(d) ap = ap' for all configurations a > e. 

Proof. It suffices to show that (a) =^ (b) and (c) =^ (d), since (b) =^ 
(c) and (d) =^ (a) trivially. 

(a) =^ (b): Suppose that ap = ap' for some particle configuration a, 
and let r G S{G/T) be a recurrent configuration. Since r is recurrent, 
there exists a particle configuration ti such that t° = t and ti > a. 
Writing Ti = a + ai for some Ui > 0, we obtain 

rp = T°p = Tip = ai{ap) 

= ai[ap ) = Tip =r^p = rp . 

(c) ^ (d): If cr > e, then writing a = e + r we have ap = T{ep) = 
T{ep') = ap'. D 

Definition. Rotor configurations p and p' are equivalent, denoted p = 
p, if the four equivalent conditions of Lemma 10 hold. 

From condition (c) of Lemma 10 it is immediate that = is an equiv- 
alence relation. We write the equivalence class of p as [p] . 

Lemma 11. If p = p' then rp = rp' for all particle configurations r. 
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Proof. If p = p' , then there exists a such that crp = ap' , which imphes 
cr(rp) = T^crp) = r^ap') = cr{Tp'), which imphes rp = rp'. D 

We say that a rotor configuration p is reachable from p' if there exists 
a particle configuration a 7^ such that p = ap'. We say that a rotor 
configuration p is recurrent if it is reachable from itself. The following 
lemma encapsulates the remaining results of [11] that we will need. 

Lemma 12. The following properties of a rotor configuration p are 
equivalent: 

(a) p is recurrent. 

(b) p is acyclic. 

(c) p = ep. 

Proof. The equivalence of (a) and (b) is Lemma 3.15 of [11]. The 
implication (b)^(c) follows from the well-definedness of the action 
of S{G/T) on acyclic rotor configurations (Lemma 9). To see that 
(c)^(a), note that ep = (ee)p = e(ep) is reachable from itself, hence 
recurrent. D 

See [11, Lemma 3.16] for several other conditions equivalent to being 
recurrent. 

Lemma 13. Each equivalence class of rotor configurations contains 
exactly one that is recurrent. The unique recurrent configuration equiv- 
alent to p is ep, where e is the identity element of S{G/T) . 

Proof. Let p be any rotor configuration. Then 

e{ep) = (e^)p = ep. 

Hence ep is recurrent by Lemma 12 and ep = p. So each equivalence 
class contains at least one recurrent configuration. 

For the reverse direction, suppose that p and p' are both recurrent 
and that p = p'. By Lemma 12 we have p = ep = ep' = p' . D 

As a consequence of Lemmas 9, 11, and 13, we have 

Corollary 14. Each addition operator Ea- acts as a permutation on 
the set R/ = of equivalence classes of rotor configurations. Thus the 
action of the sandpile monoid Q° on rotor configurations projects to an 
action 

S{G/T) X R/= -^ R/= 

of the sandpile group S{G jT) on equivalence classes of rotor configu- 
rations. 
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3.2. Cycle pushing. Lemma 13 gives one way to compute tlie unique 
recurrent rotor configuration equivalent to p: first compute the identity 
element e of the sandpile group of G/T, then add e{v) particles at each 
vertex f G Vq and stabilize. Note however that this is rather inefficient; 
for instance, in the case where p is already acyclic, a smart algorithm 
would recognize this fact and simply output p directly. We now describe 
a more efficient way to compute ep. The idea is to convert p into an 
acyclic configuration by successively removing cycles in the rotors. We 
call this process complete cycle pushing. 




Rotor 
Progression 



Rotor 
Regression 



Figure 1. Example of complete cycle pushing, start- 
ing from an arbitrary rotor configuration (upper left) to 
obtain an acyclic rotor configuration (lower left). Here 
the target set is T = {v^}. At each step, each rotor 
participating in a cycle (drawn in red) is regressed coun- 
terclockwise, until there are no more cycles. 

Cycle pushing is a key idea in Wilson's work on random stacks [19] 
and in recent work on fast simulation of rotor-routing [10]. Suppose the 
rotor configuration p contains a cycle C with vertices Wq, f i, f 2, . . . , fr = 
vq] that is, for all < j < r — 1, the arc pivj) points from Vj to fj+i. 

The rotor configuration Cp obtained by pushing C is given by 



Cp{v) 




iiveC, 
otherwise. 
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In other words, for each j = 0, 1, . . . , r — 1 the rotor p{vj) is regressed, 
and the other rotors remain unchanged. 

Suppose we have a sequence of rotor configurations po, pi, . . . , pm 
where for each i < m the configuration pj+i is obtained from pi by 
pushing a cycle Cj in pi, and suppose moreover that pm is acychc. We 
say that pm is obtained from po by complete cycle pushing. 

Complete cycle pushing involves a choice of ordering in which to 
push the cycles Cj. Wilson ["'^j showed that these choices do not affect 
the outcome: if p, p' are acyclic configurations that can be obtained 
from Po by complete cycle pushing, then p = p'. We will not use the 
uniqueness in our proofs: in fact. Lemma 17 below gives another proof 
of Wilson's result. 

Figure 1 shows how a rotor configuration is affected by cycle pushing. 
The shaded vertex f 5 is the target vertex, and the cycles that are pushed 
(first the 3-cycle t>i — )■ t>3 — > f4 — > fi, then the 2-cycle Wi — )■ ^2 — )■ fi, 
and then the 2-cycle V2 -^ v^ ^ V2) are shown in red. 

Lemma 15. If p' is obtained from p by cycle pushing, then p' = p. 

Proof. Let p' = Cp, and let a = Ic be the particle configuration con- 
sisting of one particle at each vertex Vj of the cycle. We claim that 
ap' = ap. Starting from p', let each particle take a single step of rotor 
walk: for each j, the particle at Vj moves to fj+i (taking indices mod 
r), and the rotor at Vj progresses to p'{vj)^ = p{vj). Since each ver- 
tex Vj on the cycle sends one particle to f j+i and receives one particle 
from fj-i, the resulting particle configuration is still cr; on the other 
hand, the rotor configuration has changed from p' to p. By the abelian 
property we conclude that ap' = ap, and hence p' = p. □ 

Lemma 16. For any initial rotor configuration, any sequence of cycle 
pushing moves yields an acyclic configuration in finitely many steps. 

Proof. Recall that target vertices do not have rotors. Hence if a vertex 
w has an arc to a target vertex t, then w can participate in only a finite 
number of cycle pushing moves, because at some point the rotor at w 
would point to t, and thereafter w cannot belong to a pushable cycle. 
Thereafter, each vertex v that has an arc to w can participate in only a 
finite number of cycle pushing moves, because at some point the rotor 
at V would point to w, and thereafter v cannot belong to a pushable 
cycle. Continuing in this fashion, and using the strong connectedness 
of G, we see that every vertex can participate in only finitely many 
cycle pushing moves. D 
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Lemma 17. Let p be a rotor configuration. Any sequence of cycle 
pushing moves that starts from p must terminate with ep, the unique 
acyclic rotor configuration equivalent to p. 

Proof. By Lemma 16, any sequence of cycle pushing moves starting 
from p must terminate in an acyclic configuration p'. By Lemma 15 
we have p' = p. Since p' is acyclic, p' is recurrent by Lemma 12, and 
hence p' = ep by Lemma 13. D 




cycle pushing ^i 



Rotor 
Progression 



Rotor 
Regression 



Figure 2. The rotor configurations pi at top left and p2 
at top right yield the same acyclic configuration p after 
complete cycle pushing, so they are equivalent by 
Lemma 18. 



The next lemma shows that equivalence of rotor configurations is the 
refiexive-symmetric-transitive closure of the relation p ~ Cp given by 
cycle pushing. 

Lemma 18. pi = p2 if and only if there exists a rotor configuration 
that is accessible from both pi and p2 by a sequence of cycle pushing 
moves. 

Proof. If pi = p2 then epi = ep2 is accessible from both pi and p2 by 
Lemma 17. Conversely, if p' is a configuration accessible from both pi 
and p2, then pi=p' = p2 by Lemma 15. D 



LOCAL-TO-GLOBAL PRINCIPLES 17 

For a pictorial example, see Figure 2. The rotor configuration p at the 
bottom is acyclic, and the other two rotor configurations lead to p after 
a single cycle pushing move (in one case, the 3-cycle vi -^ v^ -^ V2 -^ f i 
is pushed, and in the other case, the 2-cycle f 4 — )• ^5 — )■ v^ is pushed). 
As in Figure 1, rotors progress by turning clockwise and regress by 
turning counterclockwise. Lemma 18 tells us that the two non-acyclic 
rotor configurations pi and p2 must be equivalent, and indeed the reader 
can check that condition (a) of Lemma 10 is satisfied if one takes a to be 
the particle configuration with a single particle at W4; that is, if we add 
a single particle at v^ and let it perform rotor walk until reaching the 
target vertex fg, then the two rotor configurations become the same. 

Denote by t^{p) the target vertex reached by a particle started at v 
if the initial rotor configuration is p. 

Lemma 19. If pi = P2, then t^(pi) = ty{p2) for all f G Vq. 

Proof. By Lemma 18 it suffices to consider the case where p2 is obtained 
from pi by pushing a cycle vo,vi, . . . ,Vr = vq. If the particle added to 
pi at V never hits the cycle, then the particle added to p2 at v will 
traverse the exact same path, arriving at the same target. On the 
other hand, suppose the particle added to pi at v hits the cycle, say 
at Vq. Then the particle added to p2 at v will take the same walk to 
Vq and will then traverse the cycle, arriving back at Vq. At this point 
the rotor configuration will be the same as the rotor configuration for 
the first situation (i.e., starting from pi) when the particle first hits Vq. 
Thereafter, the two processes evolve identically, since in both situations 
the particle is at vq and the rotor configurations at this stage are the 
same in both evolutions. In particular, the particle will end up at the 
same target vertex. D 

3.3. Proof of the periodicity theorem. We can now prove our first 
main result, that the hitting sequence associated with a (periodic) rotor 
mechanism is periodic. 

Proof of Theorem 1. Let ti,t2, . . . be the hitting sequence for initial 
rotor configuration po, and for n > 1 let p„ be the rotor configuration 
after n particles released from the source vertex s have hit the targets 
ti, . . . ,tn (staying put after each hit). Then p„ = EgPn-i for n > 1. 
Let [pn] denote the equivalence class of p„. Recall that Eg acts as 
a permutation on equivalence classes (Corollary 14), so the sequence 
[po], [pi], [P2], ... is periodic, say with period D. Then by Lemma 19, 
since [pn-i-z)] = [Pn] for all n > 0, we conclude that tn+o = tn for all n > 
1, which shows that the hitting sequence is periodic with period D. D 
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Next we identify the period D of the sequence [po], [pi], . . . arising 
in the proof of Theorem 1. This in turn gives an upper bound on the 
period of the hitting sequence ti,t2,. . ., namely, the latter period is a 
divisor of D. Denote by Sg the particle configuration consisting of 1 
particle at the source vertex, and let gs = {6s + e)° be the corresponding 
recurrent configuration. 

Lemma 20. Let D he the order of Qs in the sandpile group S{G/T). 
The sequence of equivalence classes of rotor configurations {[p„]}„>o 
has period D. Moreover, the hitting sequence satisfies tn+n = tn for all 
n>l. 

Proof. For any rotor configuration p, since p = epwe have by Lemma 1 1 

gsp = (Ss + e)°p = {5s + e)p = Es{ep) = E^p. 

Since g^ = e, we obtain 

E^P = g^P = ep = p 

which shows that Pn+n = Pn for all n > 0. Conversely, if p„_|_fc = p„ for 
some n > and k > 1, then g^pn = Pn, which implies g^ = e since the 
action of S{G/T) on equivalence classes of rotor configurations is free; 
hence k must be divisible by D. 

The fact that tn+o = tn for all n > 1 follows from Lemma 19. D 

4. Time reversal and antiparticles 

4.1. Stack flipping. Recall the stacks picture introduced in §2. Each 
vertex v & Vq has a stack p^, which is a bi- infinite sequence of arcs 

Pv = [■■■, e_2, e_i, eo|ei, 62, 63, . . . ]. 

(We abuse notation slightly by using the same letter (p) to denote a 
stack configuration and its corresponding rotor configuration.) The 
Cj with i < constitute the "past" of the stack, the Cj with i > 
constitute the "future" of the stack, cq is the retrospective state of the 
stack, and ei is the prospective state of the stack; the pointer "|" marks 
the divide between past and future. When a particle at v takes a step, 
the pointer shifts to the right, so that the stack at v becomes 

[. ..,e_2,e_i,eo,ei|e2,e3,...] 

and the particle travels along arc ei. 

Shifting the pointer at v to the right corresponds to progressing the 
rotor at f , or in stack language, popping the stack at v; correspondingly, 
shifting the pointer at v to the left will be called regressing the rotor 
or pushing the stack at v. When we perform cycle pushing, the pointer 
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for the vertex v moves one place to the left for all vertices v belonging 
to the cycle. 

We define stack flipping as the operation on a bi-infinite stack that 
exchanges past and future, turning 

[. ..,e_2,e_i,eo|ei,e2,e3,...] 

into 

[. ..,63,62, ei|eo,e_i,e_2,...]. 

Given a stack configuration p = (pj,)^gvb, let $(p) denote the stack 
configuration obtained by flipping all its stacks. Note that <l>(<l>(p)) = p. 

Lemma 21. Let p be a rotor configuration that has a cycle C. Then C 
is also a cycle of^{Cp), and 

Hcmcp))) = p. 

Proof. Let f be a vertex of C. Let p' = Cp, and write the rotor stack 
at V as 

p^ = [..., 6_2, 6_i, 6o|6i, 62, 63, . . . ]. 

If we push the cycle, the stack at v becomes 

(Cp)„ = [. . . ,6_2,6_i|6o,6i,62,63,. ..]. 

If we then flip all the stacks, we obtain 

{^(Cp))^ = [..., 63, 62, 61, 6o|6_i, 6_2, . . . ]. 

The retrospective rotors at the vertices f G C are now as they were 
initially in p, so they form the same cycle C. Pushing that cycle yields 

(C(<I>(Cp)))^ = [. ..,63,62, 6i|6o,6_i,6_2,...]. 

Finally, flipping the stacks once more brings us to 

(<I>(C(<I>(Cp))))^ = [•••, e_2, e_i, 6o|6i, 62, 63, ... ] 

which equals p„. 

Meanwhile, for those vertices v that are not part of the cycle C, the 
stack at v is simply reversed twice (with no intervening cycle pushing 
moves to complicate things), so this stack ends up in exactly the same 
conflguration as it was in p. D 
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Diagraminatically, writing p' = Cp, we have 

c 

P ^ P 



$ 



$ 



Hp) —^ Hp) 

Note the reversal of the direction of the C arrow. 

Lemma 22. If p = p' , then $(p) = $(p')- 

Proof. By Lemma 18, it suffices to show that if two stack configurations 
are related by a cycle pushing move, then their fiips are related by a 
cycle pushing move. But that is precisely what Lemma 21 tells us. D 

4.2. Antiparticles. Next we introduce antiparticles. Like particles, 
they move from vertex to vertex in the graph, but they interact with the 
stacks in a different way. Suppose that the current stack configuration 
at V is 

[. ..,e_2,e_i,eo|ei,e2,e3,...] 

and that there is an antiparticle at v. An antiparticle step consists of 
first moving the particle along the arc Cq and then pushing the stack 
at V to obtain 

[...,e_2,e_i|eo, 61,62, 63,...]. 

(Compare: a particle step consists of first popping the stack at v to 
obtain 

[. ..,6_2,6_i,6o,6i|62,63,...] 

and then moving the particle along the arc 6i.) 

Lemma 23. If p' is obtained from p by moving a particle from v along 
arc e, then $(p') is obtained from $(p) by moving an antiparticle from 
V along arc e. 

Proof. Write the stack at v for the rotor configuration p as 

p^ = [..., 6_2, 6_i, 6o|6i, 62, 63, . . . ]. 

When a particle at v advances by one step, the particle moves along 
the arc 61 and the stack at v becomes 

P'„ = [• • • , e-2, 6_i, 60, 6i|62, 63, . . . ]. 

On the other hand, the stack at v for the fiipped rotor configuration 
$(p) is 

^{P)v = [•••,63,62,6i|6o,6_i,6_2,...]. 
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When an antiparticle at v advances by one step, the antiparticle moves 
along the arc ei and the stack at v becomes 



.e3,e2|ei,eo,e_i,e_2, 



which equals $(p')«- 



D 



Just as one defines particle addition operators E^, one can define 
antiparticle addition operators E^ on rotor configurations: to apply 
E~ , add an antiparticle at v and let it perform rotor walk on G (using 
the antiparticle dynamics described above) until it arrives at a vertex 
in the target set T. To highlight the symmetry between particles and 
antiparticles we will sometimes write E^ instead of E^ for particle 
addition operators. Note that in general, E^ and E~ do not commute. 

Write tj(p) (resp. t~{p)) for the target vertex hit by a particle (resp. 
antiparticle) started at v if the initial rotor configuration is p. 

Lemma 24. For any rotor configuration p and any v & Vq we have 
HE^ip)) = E~mp)), and t+(p) = t-($(p)). 

Proof. This follows by repeated application of Lemma 23: the sequence 
of vertices traveled by the particle added to p at u is the same as the 
sequence of vertices traveled by the antiparticle added to $(p) at v. D 

Diagrammatically, writing p' = E^p, we have: 

E+ 



P 



-* P 



$ 



<|.(p) 



E- 



$ 



$(p') 



Lemma 25. If p = p' , then E^ p = E^ p' and t^ (p) = t^ (p') for all 
veVo. 

Proof. We have E~ = $ o E^ o $ by Lemma 24. Moreover, $ preserves 
equivalence by Lemma 22 and E^ preserves equivalence by Lemma 11, 
so E~ must preserve equivalence. This proves the first statement. For 
the second, since $(p) = ^{p'), we have by Lemmas 24 and 19 



t.-(p) = ttmp)) = ttmp')) = t-{p'). D 
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4.3. Loop-erasure. If a path (xq, ... ,a;,.) in the directed graph G con- 
tains a cycle, i.e., a sub-path (xp, Xp+i, . . . , Xq) with Xq = Xp, define the 
first cycle as the unique cycle with q as small as possible; we may re- 
place the path by the shorter path (xq, . . . , Xp-i,Xp, x^+i, . . . , Xr) from 
which the q — p vertices of the first cycle have been removed. If this 
new path contains a cycle, we may erase the first cycle of the new path, 
obtaining an even shorter path. If we continue in this fashion, we even- 
tually obtain a simple path from xq to Xr, called the loop- erasure of 
the original path. 

The notion of loop-erasure is due to Lawler [15], who studied the 
loop-erasure of random walk. As is mentioned at the end of §5 of [ ], 
there is also a connection between loop-erasure and rotor walk. Given 
a rotor configuration p and a set 5* C Vq, define popping S as the 
operation of popping the stack at each vertex in S to obtain the new 
rotor configuration 



S+p{i 




iiv eS, 

otherwise. 



(Compare to cycle pushing §3.2, in which the rotors are regressed in- 
stead of progressed.) For a rotor configuration p and a vertex v E Vq, 
let r be the loop-erasure of the path xq, . . . ,Xr traveled by a particle 
performing rotor walk starting from xq = v until it hits the target set. 
Let Ci, . . . , Cm be the cycles erased to obtain F. For any vertex w, the 
number of z (0 < i < r — 1) with Xi = w is equal to the number of j 
(1 < j < m) for which w E Cj, plus either 1 or according to whether 
or not w &r. Hence the final rotor configuration E^p can be obtained 
from p by popping the cycles Ci, . . . , Cm and the path 7 := F — {x^}; 
that is, 

E> = 7+C+...0- (1) 

Lemma 26. For every rotor configuration p and every v E Vq we 
have E^ E^ p = p, and the path traversed by the antiparticle is the 
loop-erasure of the path traversed by the particle. In particular, the 
antiparticle hits the same target as the particle: 

t-{Etp)=tt{p). 

Proof. After the particle has been added to p at v^ changing the rotor 
configuration to E^ p and arriving at target t = t^{p), the retrospective 
rotor at each vertex v is the arc that the particle traversed the last time 
it left V. Hence the rotors of E^p give a simple (cycle-free) path 7 from 
V to t, and the antiparticle will travel this path, arriving at the same 
target t. By (1), the rotor configuration E^p is obtained from p by a 
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sequence of cycle-popping moves followed by a "path-popping move" 
along 7. The motion of the antiparticle from v to t undoes the path- 
popping move, so all that survives in E~ E^ p are the cycle popping 
moves. Since cycle popping doesn't change the equivalence class of a 
rotor configuration (by Lemma 18), we conclude that E~E:^p = p. D 

Likewise, for every p we have E^E~ p = p. Lemma 26 thus says 
that the products E:^E~ and E~E^ act as the identity operation on 
equivalence classes of rotor configurations. That is, if we view E^ and 
E~ as elements of the sandpile group (acting on equivalence classes of 
rotor configurations), they are inverses. 

4.4. Proof of the rotor-reversal theorem. Now we turn to the 
proof of our second main result, that reversal of the periodic pattern 
of the rotor mechanism at all vertices causes reversal of the periodic 
pattern of the hitting sequence. To save unnecessary notation in the 
proof, we write E^ := Ef and t^ '■= tf. 

Proof of Theorem 2. As in the proof of Theorem 1 , the sequence of 
equivalence classes, [po], [pi], [P2], • • • is periodic, say with period D. 
Now consider the hitting sequence for antiparticles released from the 
source vertex s from initial configuration po- Define 770 = po and rji = 
E^ijii^i) for z > 1. We first show by induction on i that r/j = po-i for 
alH = 0, . . . , D. The base case i = is the fact that po = Pd] and for 
1 < z < D, if rji^i = pD-i+i then by Lemmas 25 and 26, 

Tji = E-{r]i_i) = E'{pD-i+i) = E~E~^{pD-i) = pD-i 

which completes the inductive step. 

Now for i > I, let u^ = t+(pj_i) and u~ = t~{r]i_i) be the hitting 
sequences for a particle (resp. antiparticle) started at s with initial rotor 
configuration po. Using the second statements of Lemmas 25 and 26, 
we have for z = 0, . . . , D — 1 

M^+i = t-{ri.i) = r{pD-i) = t-{E+pD-i-l) = t-^{pD-t-i) = u^_^. 

By Lemma 24, the hitting sequence for a particle starting at s with 
rotor configuration $(po) equals the hitting sequence for an antiparticle 
starting at s with rotor configuration po, that is, the sequence {Mj~}j>i. 
Moreover, since rjD = pq = pd = tjq, the sequence {u~}i>i satisfies 
u~^^ = u~ for all i > 1 by Lemma 25. Hence the particle hitting 
sequences for $(po) and po are both periodic modulo D, and reversing 
the first D terms of the latter hitting sequence yields the first D terms 
of the former. D 
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